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Nuclear condensation and symmetry energy of dilute nuclear matter: 

an S-matrix approach 

J.N. DcQ and S.K. SamaddaiQ 
Theory Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata 700064- 

Based on the general analysis of the grand canonical partition function in the S-matrix framework, 
the calculated results on symmetry energy, free energy and entropy of dilute warm nuclear matter 
are presented. At a given temperature and density, the symmetry energy or symmetry free energy of 
the clusterized nuclear matter in the S'-matrix formulation deviates, particularly at low temperature 
and relatively higher density, in a subtle way, from the linear dependence on the square of the 
isospin asymmetry parameter X = (p„ — Pv)l {Pn + Pp) j contrary to those obtained for homogeneous 
nucleonic matter. The symmetry coefficients, in conventional definition, can then be even negative. 
The symmetry entropy similarly shows a very different behavior. 

PACS numbers: 12.40.Ee, 21.65.Mn, 21.65.Ef, 24.10.Pa 

Keywords: symmetry energy, symmetry entropy, nuclear matter, statistical meehanies, S-matrix 



I. INTRODUCTION 

Understanding properties of symmetry energy for low 
density nuclear matter is of much topical interest in both 
astrophysical and laboratory context. The supernova 
simulation dynamics has a sensitive dependence on the 
symmetry energy a higher symmetry energy, for ex- 
ample, leads to a lower electron (e~)-capture rate in the 
supernova collapse phase that may result in a stronger 
explosive shock. The variations in the e~-capture rate 
also produces changes in the neutrino luminosities that 
are potentially observable. The isotopic abundance of 
relatively heavier elements in explosive nucleosynthesis is 
further directly correlated to the symmetry energy. The 
neutron skin thickness of heavy nuclei has also a direct 
dependence on the density variation of the symmetry en- 
ergy [1,3. At subnormal densities (0.2 < p/po < 1.0, 
Po is the saturation density of nuclear matter) , the sym- 
metry coefficients that are a measure of the symmetry 
energy per baryon {esym) have recently been estimated 
from analysis of data related to isotopic distributions [J , 
isospin diffusion Q , and isoscaling 0, Q in heavy ion 
reactions. The different analyzes give somewhat differ- 
ent results. In the density domain mentioned, the sym- 
metry coefficient Ce {=€.syml X'^) varies with density as 
Ce(p) ^ Ce{pq){pI Pa)'^ with 7 lying in a broad range 
Q, between 0.55 — 1.05. A more definitive answer 
to the question of the density dependence of symmetry 
energy is still wanting. 

On the theoretical side, symmetry energy of infinite 
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nuclear matter has been calculated over a wide density 
range in many guises; many-body theories using various 
nucleon-nucleon interactions or interaction Lagrangians 
have lead to varying results [l3|- There have also 
been some recent investigations on the symmetry free 
energy of hot nuclear matter [ll| and also of finite nuclei 
p^ . These calculations differ in details, but the general 
qualitative behavior of the symmetry coefficients with 
density do not deviate much from the experimental trend. 
A similar power law variation is exhibited; the exponent 
7 lies mostly within the broad limits as extracted from 
experimental analyzes. 

The above calculations have been done in the mean- 
field (MF) approximation; the system is taken to be ho- 
mogeneous nucleonic matter. For dilute nuclear matter, 
however, the system minimizes its total energy or free en- 
ergy by forming clusters [l^, [ij] . A detailed knowledge 
of the composition of nuclear matter is then needed to 
appreciate how the symmetry energies are affected when 
matter gets clusterized. This has a direct role in a bet- 
ter understanding of neutrino-driven energy transfer in 
inhomogeneous supernova matter p^ . Using the virial 
expansion technique, clusterization in dilute nuclear mat- 
ter and its import in the evaluation of the symmetry 
coefficients has been investigated recently by Horowitz 
and Schwcnk where they considered the matter to 
be composed of n,p, a. The investigations with the same 
cluster species have been followed further |l3| to connect 
experimentally the temperature and density-dependent 
entropic contributions to the symmetry free energy co- 
efficient. The resulting symmetry coefficients are found 
to be considerably larger than the corresponding ones 
obtained from MF calculations. This is an important 
result, it shows the strong role of clusterization and nat- 
urally calls for a realization of symmetry coefficients if 
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all possible permissible clusters are incorporated in the 
calculation. The present article is an attempt in this di- 
rection. 

II. ELEMENTS OF THEORY 

The grand canonical partition function in the iS-matrix 
formalism of statistical mechanics proposed by Dashen et 
al sets the logical framework of our calculation. The 
details of the formalism, as applied to nuclear matter 
is given in Ref. [l^; for the sake of completeness, the 
essentials are presented in the following. 

The partition function Z for the two-component nu- 
clear matter composed of the elementary species neutrons 
and protons is written as 

Z = Tre"'^'-^"^"*""''"*"^ , (1) 

where (3 is the inverse of temperature T of the system, 
H the total Hamiltonian, Np_n the number operators for 
protons and neutrons, and /Xp^„ are the corresponding 
chemical potentials. The partition function can be de- 
composed as 

00 

2= E e,CTrz,Ne-^", (2) 

ZM=Q 

where the fugacities are given by Cp = e''^" , Cn = ■ 
The trace Trz,N is taken over states of Z protons and N 
neutrons. For small C,p and the quantity ln2^ can be 
expanded in a virial series 

luZ^Y^DzMpCn ■ (3) 
Z,N 

The prime on E indicates that the term with Z = A^=0 is 
excluded. Evaluation of the virial coefficients Dz.n gives 
the partition function and thence the thermodynamic be- 
havior of the system. 

Following the temperature-Green's-function method, 
it was shown in Ref. [l8| that all the dynamical informa- 
tion concerning the microscopic interaction in the grand 
potential of the system can be collected in two types of 
terms so that the partition function is written as 

lnZ = lnZ(l+lnZ,,,t. (4) 

The first term corresponds to contributions from stable 
single particle states of clusters of different sizes (neu- 
trons and protons included) formed in the infinite system 
behaving like an ideal quantum gas [the superscript (0) 
indicates this behavior] and the second term corresponds 
to contribution from multiparticle scattering states, re- 
spectively. The particle piece can further be split in con- 
tributions from ground states and excited states of the 
bound nucleon clusters, so that 

lnZ^l=lnZ(0)+lnZ(0). (5) 



The first term in Eq. (5) is a sum of ideal gas terms, one 
for each of the ground states of all the possible species 
of mass A with Z protons and N neutrons that can be 
formed in the system, 

Z,N •' ^ ' 

ln(lTCz.7ve-''(^'/^^™)). (6) 

Here m is the nucleon mass, p is the momentum of the 
nucleus, and Cz.w the 'effective fugacity' given by Cz.w = 
g/3(Mz,N-i--Bz,N)_ Bz.N is the binding energy of the nucleus 
and [iz,N its chemical potential. From the condition of 
chemical equilibrium among the different species, [iz.N = 
Z[ip -I- 7V/^„. 

The =F sign in Eq. (6) corresponds to nuclei with A even 
or odd, obeying Bose or Fermi statistics, V is the volume 
of the system and g is the spin degeneracy. The sum 
includes the original elementary species, namely the neu- 
trons and protons. The Coulomb interaction is assumed 
absent in nuclear matter, ideally the sum in Eq. (6) would 
involve infinite terms as the maximum cluster mass A can 
then even be infinite. However, for applications to real 
physical systems such as neutron star matter. Coulomb 
effect is to be included in the binding energies of nuclei. 
In that case, the sum is finite, conditioned by the sta- 
bility of nuclei with inclusion of Coulomb in the binding 
energies. Calculations in the ^-matrix formalism with 
inclusion of Coulomb in the fragment binding energies 
would henceforth be referred to as SMF, those without 
Coulomb would be called SNC. Eq. (6) can be readily 
expanded in a virial series as 

'»4? = ''E^(&»±i#+-). m 

in powers of effective fugacities provided |Cz,Af| < 1- 
Here, A(Am) = ■\/27r/ (AmT) is the thermal wavelength 
of a cluster of mass Am. We work in natural units 
h = c = 1. A nucleus in a particular excited state is 
taken as a distinctly different species and can be treated 
in the same footing as the ground state. The density of 
states is quite high in relatively heavy nuclei, increases 
nearly exponentially with the square root of excitation 
energy E and so the contribution of the excited states of 
a single nucleus is written as an integral over E of the 
ideal gas term weighted with the level density uj{A, E), 

\nZi"J = 9 '^^ '^(^' ^) ^ 

Z,N "^^0 

j^^H'^Cz^^^-''^''^'---^). (8) 
For the level density, we take the expression [2^ 

-(Ai^) = ^^. (9) 
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The level density parameter a is taken as A/8 MeV""'^, 
its empirical value. 

In Eq. (8), the prime on E denotes exclusion of the 
light nuclei {A < 8) from the sum. For these nuclei, 
we take only the ground states; their degeneracy factor 
g is taken from experiments. For other nuclei, for both 
ground and excited states, g is taken as 1 or 2, according 
as they are bosonic or fermionic. The lower limit Eq is 
dictated by the location of the first excited state. The 
upper limit Es is the separation energy. We take Eq=2 
MeV and Es =8 MeV. 

The scattering piece In Zgcat of Eq. (4) can be formally, 
but explicitly written (loj for nuclear matter: 



In Zscat= V 



E 



Here, the double sum refers to the sum over all possible 
scattering channels, each having its chemical potential /i 
and formed by taking any number of particles from any 
of the stable species (proton, neutron and nuclei in their 
ground and excited states) and the trace is over all plane 
wave states for each of these channels. S is the scattering 
operator and A the boson symmetrization or fermion an- 
tisymmetrization operator. The subscript c denotes only 
the connected parts in the diagrammatics of the expres- 
sion in the parenthesis. A channel in the set has a total 
number of Zt protons and Nt neutrons {At ~ Nt + Zt); a 
denotes all other labels required to fix a channel within 
this set. Bzt.Nt,a- is the sum of the individual binding 
energies of all the particles in the channel and e is the 
total kinetic energy in the cm. frame of the scattering 
partners. Examination of Eq. (10) shows that channels 
with larger binding energies are more important, because 
of the factor e'^^^t "*-" . Furthermore, two-particle chan- 
nels are expected to dominate over multiparticle channels 
with the same Zt and Nt from binding energy consider- 
ation. We therefore consider only two-particle scattering 
channels. It is convenient to divide the channels into light 
ones, consisting of low mass particles {A < 8, say) and 
heavy ones {A> 8), so that we write 



In Z., 



InZ, 



light 



In 2, 



heavy ; 



(11) 



as the sum of contributions from the light and heavy 
channels. 

Experimentally it is known that the scattering of rela- 
tively heavier nuclei is dominated by a multitude of res- 
onances near the threshold. The S'-matrix elements are 
then approximated by resonances. Following [2l|, [2^ . 
each of these resonances are treated like an ideal gas 
term and then In Zheavy can be written in the same form 
of ^ assuming their level density to be the same 

as those of the excited states. The integration over E in 
Eq. (8) now extends from Eg to i?,., the limit of resonance 
domination. The damping of the integral in Eq. (10) due 



to the presence of the Boltzmann factor assures contri- 
bution only from low energies; we therefore take iJ^ — 12 
MeV. 

The scattering channels NN, Nt, NHe^, Na {N refers 
to nucleon) and aa are considered for evaluation of 
Zught- Inclusion of other light particle scattering chan- 
nels may have some influence on the present results. In 
intermediate energy heavy ion collisions, for example, 
deuteron usually has considerable multiplicity and the 
contribution to In Zught from deuteron scattering chan- 
nels is worth further exploration. With the choice of the 
scattering channels as mentioned above, we then write 



InZ, 



light 



InZ 



NN- 



InZ^rt-l-lnZ 



NHe^ ' 



AnZ 



No 



~ln Zn 



(12) 

We explicitly write the contribution from the NN chan- 
nel. If we consider only elastic two-body scattering, the 
trace in Eq. (10) becomes a sum over the derivative of 
the phase shifts of the appropriate partial waves. It gives 
formulas of the same form as derived by Beth and Uh- 
Icnbeck [2^ for the second virial coefficient. It is given 
as 



In Znn = 



V 



A3(2m) 

+CpCn(-3 + A^^i + A^l,")}. (13) 



Here / is the isospin index and 
1 1"°° 

KN^-f dee-^^Y.(^J+^)^2slir (14) 

S,L,J 

The quantity (e) refers to the NN phase shift in 

the LSJ channel. The contributing partial waves are de- 
termined by / through the requirement of antisymmetry 
on the total wave function of the NN system. The other 
terms in Eq. (12) have nearly similar forms [l^. The A's 
for the NN, Na and aa channels are evaluated in [l^ 
and those for NHe^ and Nt are available in [23|. This 
completes the evaluation of partition function for nuclear 
matter. It is then straightforward to get the relevant ob- 
servables like the pressure P, the number density of the 
i-th species pi, free energy per baryon / or the entropy 
per baryon s from the relations. 



P = T- 



,lnZ 



Pi = Q 



/ d InZ 
\dQ~Jv 



with the baryon density p — AiPi. 

III. RESULTS AND DISCUSSIONS 

We have calculated the nuclear equation of state 
(EOS), fragment distributions, the symmetry free en- 
ergy, symmetry entropy and the symmetry coefficients 



(symmetry energy coefficient C e and the symmetry free 
energy coefficient Cp) for dilute nuclear matter. The 
calculations have been done in the SMF and the SNC 
model where all possible nuclear species are considered 
and compared, in order to explore the role of heavy 
species, with calculations [l^ [2^ where only the light 
species (n,p^ d, t, He^ and a) are taken. The latter model 
is subsequently referred to as the light species model 
(LSM). The three models (SMF, SNC and LSM) would 
be collectively referred to as the condensation models. To 
highlight the importance of clusterization on the physical 
observables, the MP results are also presented. 

In practice, an asymptotic wave function may not have 
a precise meaning at relatively high density; it would then 
be difficult to have a meaningful expression of the parti- 
tion function in terms of the S'-matrix elements. We, 
therefore, restrict our calculations to low density nu- 
clear matter and have considered up to a baryon density 
p=0.01 fm-3. 

In Fig. 1, the EOS (P — p) is displayed for symmetric 
nuclear matter (p„ ~ pp) at temperatures T =2, 4, and 
8 MeV. The dotted lines refer to calculations in the MF 
model with the SkM* interaction. At low temperatures, 
it is seen that the system enters the unphysical region 
{dP/dp < 0) in the density range considered. This can be 
avoided by applying Maxwell's construction. In the SMF, 
because of the many-body correlations (condensation), 
this unphysical behavior does not arise. Since Coulomb 
interaction is absent in nuclear matter, to compare re- 
sults from the mean- field, a set of S'-matrix calculations 
have been done with Coulomb switched off in the nuclear 
binding energies (SNC). The SNC calculations are repre- 
sented in the figure as dot-dashed lines. With isothermal 
compression, at lower temperatures, the pressure levels 
off at very low densities as shown by the dot-dash lines in 
Fig .1(a) and 1(b), signaling a behavior like a first-order 
phase transition. At the higher temperature 8 MeV, the 
said transition starts at a density beyond 0.01 fm~^, it is 
not seen in the figure. The sum in Eq. (3) for the SNC 
calculation runs up to infinity in principle; in practice, 
one takes a finite sum for calculational facilitation. We 
have taken the maximum mass Amax=WOO. The results 
are found to be not very sensitive to further increase of 
Amax [11]. The binding energies of these nuclei are ob- 
tained using a simple liquid-drop mass formula (26j with 
Coulomb switched off. 

To help comparison with physical systems such as neu- 
tron star matter, we have also considered phenomenologi- 
cal binding energies (Coulomb included) [23] of the nuclei 
that limits the number of terms in the sum to ~9000 nu- 
clei in their ground states with Amax= 339 and Z^ax 
=136. The EOS in the SMF are shown by the full fines. 
It is seen that with Coulomb in the binding energies, the 
signature of the first-order phase transition is washed out 
with monotonic increase in pressure on isothermal com- 
pression. The results in the LSM model are shown by 
the dashed lines. At a given density, in the SMF, the 
fragment multiplicity is comparatively lesser as more 
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nucleons get bound in larger clusters; the pressure is also 
then lesser compared to the LSM model; the deviations 
are significant, particularly at higher density and at lower 
temperature. 

The composition of matter at different densities and 
temperatures for symmetric nuclear matter are shown in 
Fig. 2 through the charge multiplicities. The left, mid- 
dle and right panels correspond to temperatures T=2, 4, 
and 8 MeV, respectively, at three different baryon den- 
sities p = 0.0001, 0.001 and 0.01 fm'^. The fuU lines 
display the results from the SMF, the dashed lines are 
the ones from the LSM model and the dot-dashed lines 
correspond to the ones from the SNC model. Examina- 
tion of the results brings out a few important findings: 
i) At very low densities, the multiplicities up to Z =2 
are practically the same in the SMF and LSM models. 
In the SMF, heavier fragments may be formed, but that 
is insignificant, ii) With increase in density, heavy frag- 
ment formation can no longer be neglected. Increase in 
temperature hinders the heavy cluster formation. The 
multiplicity distributions display a saw-toothed nature. 
This is odd-even effect due to inclusion of pairing in the 
phenomenological binding energies of the nuclear clus- 
ters; this effect is diluted with increase in temperature, 
iii) At relatively low temperatures and higher densities, 
in the SNC model, the matter consists of only nucleons 
and very heavy nuclei; the matter resembles liquid-like 
along with a negligible fraction of nucleonic gas. These 
features are observed at all the three densities at T =2 
MeV and at p =0.01 fm^^ ^j. y ^4 y^^y^ These results 
are not shown in the figure. With increasing temperature 
and decreasing density, the liquid-like structures disap- 
pear. These features appear from a delicate dependence 
of the chemical potential on the density and temperature 
and the dependence of fragment binding energy on in- 
creasing fragment mass which tends to saturate at ^16 
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FIG. 1: (Color online) The EOS for symmetric nuclear matter 
at r= 2, 4, and 8 MeV. Calculations are shown for models 
in mean-field (dotted cyan), LSM (dashed black), SNC (dot- 
dashed magenta), and SMF (full line, blue) 




FIG. 2; The charge distributions at different temperatures 
and densities as shown are compared in the models of LSM 
(dashed hne), SNC (dot-dashed) and SMF (full line). 



FIG. 4: (Color online) The symmetry entropy Ssym as a 
function of X'^ at different temperatures and densities. The 
notations are the same as described in the caption to Fig. 3. 



MeV per nucleon. 

The evolution of symmetry free energy per baryon fsym 
as a function of asymmetry X'^ at different temperatures 
and densities is displayed in Fig. 3. The symmetry free 
energy for a given density and temperature is defined as 

Um = f{X)-fiO) + -Y,[P^{X)~P^mBl, (16) 
^ i 

where Bl is the Coulomb contribution to the binding en- 
ergy of the i-th fragment species. The symmetry energy 
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FIG. 3: ( Color on-line) The symmetry free energy shown as 
a function of at different temperatures. The dashed, dot- 
dashed and full lines correspond to p =0.0001, 0.001 and 0.01 
fm^"^, respectively. The results are calculated in the models 
of the mean-field (thin cyan), LSM (thick black), SNC (thick 
magenta), and the SMF (thin blue). 



Gsym can be defined Ukewise. The dashed, dot-dashed 
and the full lines correspond to calculations at densi- 
ties of p =0.0001, 0.001 and 0.01 fm-^, respectively. In 
Figs. 3(a), 3(b) and 3(c), the symmetry free energies per 
nucleon in the LSM (thick black lines) and the SMF (thin 
blue lines) arc compared at the three densities at differ- 
ent temperatures. As seen in Fig. 3(a) at T =2 McV, at 
the lowest density, the two calculations yield nearly the 
same results; with increasing density, the difference in the 
two model predictions shows up prominently because of 
the formation of larger clusters in the SMF. This differ- 
ence washes out gradually with increasing temperature 
and the results are not discernible as seen in Fig. 3(b) 
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FIG. 5: (Color online) The symmetry coefficients Ce and 
Cf as a function of density at different temperatures in the 
models of mean-field (dotted cyan), LSM (dashed black), SNC 
(dot-dashed magenta) and SMF (full lines, blue). 
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FIG. 6: The symmetry coefRcients Ce and Cf displayed as 
a function of temperature at densities p=0.0001, 0.001, and 
0.01 fm""^. The notations are the same as described in the 
caption to Fig. 5. 



and Fig. 3(c). In Fig. 3(d), the symmetry free ener- 
gies at T =2 MeV for the three densities from the SNC 
(thick magenta Unes) and MF (thin cyan hnes) models 
are compared. In the liquid-drop mass formula, the sym- 
metry energy is taken to be linear in X^. This is seen to 
be nearly true also for symmetry free energy of nuclear 
matter in a density region around the saturation density 
[m in the MF model; we find the same for dilute nuclear 
matter as is shown in Fig. 3(d). The symmetry energies 
can then be written as 



fsym — Cp X"^ 



(17) 



where Ce and Cp are the symmetry energy and sym- 
metry free energy coefficients. Clusterization affects this 
linearity at low temperatures; as the temperature is in- 
creased, the linearity tends to be restored as seen from 
Fig. 3(b) and Fig. 3(c). The symmetry energy has a sim- 
ilar behavior and is not shown here. 

One interesting result that is borne out of the SMF cal- 
culation is that at a given T and p, asymmetric nuclear 
matter may become more stable than symmetric matter 
[f{X) < /(O)]. This may result in negative fsym as is 
shown in Fig. 3(a). This happens at relatively higher 
densities. We find it to be due to the absence of isospin- 
conjugate (mirror) nuclei because of Coulomb interac- 
tion. Restricting the sum in Eq. (3) to only mirror nuclei 
removes this negativity. In the SNC model, clusters al- 
ways occur in isospin-conjugate pairs, fsym is then always 
positive as seen from Fig. 3(d). 

The symmetry entropy per baryon, defined as 
Ssym = s{X) — s(0) is presented as a function of asym- 
metry for different densities and temperatures in Fig. 4. 
The notations used are the same as those used for Fig. 3. 
Comparison between the MF and the SNC models at the 
three densities is displayed in Fig. 4(b) and Fig. 4(d). 



In the MF model, the symmetry entropy decreases with 
asymmetry. This is akin to the entropy of mixing . The 
entropy per nucleon of an ideal two-component nucleon 
gas is given, from Gibbs-Duhem relation, by 



six) = ^_^lnC„-^lnCp. 
2 p p 

Using the fact that for low density nucleonic matter 



(18) 



Pn,p — ^3 Cri.pj 



(19) 



the symmetry entropy can be shown to behave as 

SsymiX) = -^X^(^l + ^X^ + ...}j , (20) 

which for low values of asymmetry decreases linearly with 
X^. This is independent of both temperature and den- 
sity and is nearly manifested for the low density nuclear 
matter we have considered. As can be seen from the com- 
parison with the SNC calculations, clusterization changes 
this behavior; in the density and temperature domain 
where clusterization becomes important, the symmetry 
entropy is larger compared to that in the MF model and 
can even be positive. In Fig. 4(a) and Fig. 4(c), com- 
parison is made between the results from LSM and SMF. 
At high temperatures and low densities, the difference in 
the results from these models is not discernible; however, 
at relatively high density and low temperature, Ssym in 
the LSM is appreciably larger than that in the SMF as 
displayed in Fig. 4(a) for p =0.01 fm'^ and T = 4 MeV. 
This is attributed to the relatively rapid growth of mul- 
tiplicity (mostly of neutrons at the cost of p, He^ and a) 
in LSM model compared to that in SMF. 

The conventional definition of the symmetry coeffi- 
cients Ce and Cp as given by Eq. (17) holds good when 
Ssym and fsym are linear or nearly linear in X^. In Fig. 5, 
we display these symmetry coefRcients as a function of 
density for temperatures when the symmetry energies are 
nearly linear in X^ in the different models. In the den- 
sity region investigated, in the MF model, the symmetry 
coefficients increase linearly with density in contrast to 
a power-law dependence ~ (p/po)'^ with 7 ^^0.7 at rela- 
tively higher densities [l^ . The symmetry coefficients are 
enhanced significantly when clusterization is considered. 
This is most prominent in the SNC model (dot-dashed 
lines), where as already stated, at lower temperatures 
and higher densities, the system is more liquid-like; this 
is reflected in the value of Ce in Fig. 5(a), which at higher 
densities saturates to ~28 MeV, the value for normal nu- 
clear matter as given in the liquid-drop mass formula we 
use [2^. The symmetry coefficient Ce in the LSM and 
in SMF are close at high temperature T= 8 MeV, but 
differ at T= 4 MeV where clusterization is more impor- 
tant. This difference is filled up for symmetry free energy 
coefficient Cp due to the symmetry free entropy [shown 
in Fig. 4(a)] resulting in practically the same Cp in both 
the models. 
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The temperature dependence of the symmetry coeffi- 
cients Ce and Cp at a few densities in aU the models 
considered are displayed in the left and right panels of 
Fig. 6, respectively. In the mean-field model, shown by 
the dotted lines, Cb is linear and a ppr oximately constant 
for a given density as seen earlier [28|. At very low den- 
sity p=0.0001 fm~^, it is nearly zero and not discernible 
in the figure. In the condensation models a pronounced 
increase in the values of Ce is observed, particularly at 
low temperatures. In the SNC model this value is close 
to that for normal nuclear matter as expected. At larger 
temperatures, values oi Ce obtained in the condensation 
models approach those calculated from MF. 

The symmetry free energy coefficient Cp, by our 
choice, can be written as 

Cf = Ce~t'-^. (21) 

For dilute nuclear matter, in the MF model, Ce is linear 
in T at a particular density and Ssym, as stated earlier, is 
negative and proportional to X^. A linear increase of Cp 
with temperature is therefore expected in the mean-field 
model which is realized as displayed by the dotted lines 
in the right panels of Fig. 6. At higher temperatures, re- 
sults for Cp in all the models tend to merge, particularly 
at low density because of the hindrance to form clusters. 
At lower temperatures, Cp is significantly higher with 
clusterization compared to that in the MF. In Eq.(21), 
the first term on the right hand side decreases with tem- 
perature whereas the second term increases with conden- 
sation; this interplay causes a minimum in Cf which is 
more pronounced at lower densities. 

IV. CONCLUDING REMARKS 

The role of condensation on some properties of dilute 
nuclear matter, namely, the nuclear EOS, the symme- 
try free energy, entropy and the symmetry coefficients 
has been addressed in this article in the 5-matrix frame- 
work. This approach has the advantage that the relevant 



observables can be directly connected to the experimen- 
tally measured quantities like the nuclear binding ener- 
gies and the scattering phase-shifts. This approach con- 
tains no interaction potential parameters and hence, the 
results are mostly model- independent. 

Except at very low densities, the nuclear EOS in the 
S'-matrix approach differs appreciably from the one ob- 
tained in the MF model. The MF model, supplemented 
with Maxwell's construction displays a liquid-gas phase 
transition on isothermal compression, in the S'-matrix 
approach, the system responds to the compression by a 
marked growth of clusters out of the dilute nucleonic gas. 
The growth is hindered with isochoric heating. 

One of the very remarkable features of the results on 
symmetry energies (both Csym and fsym) in SMF is that 
the symmetry energies are generally nonlinear in X^, con- 
trary to the results in the MF model. The symmetry 
energies may even be negative in the SMF at low tem- 
peratures, high densities and at small X. The symmetry 
entropy similarly displays a subtle behavior with cluster- 
ization. In the conventional definition, at a particular 
temperature and density, the symmetry coefficients are 
then no longer independent of the asymmetry parameter 
and may be negative. In the region where the symme- 
try energies arc practically linear in X^, the symmetry 
energy and free energy coefficients in the S'-matrix ap- 
proach are found to be appreciably larger compared to 
those obtained in the MF model. The nuclear EOS and 
the symmetry energies and coefficients are thus seen to 
have a significant dependence on the many-nucleon cor- 
relations or cluster formation in the low density nuclear 
matter. Inclusion of these effects arc thus warranted for 
the study of physical phenomena like supernova dynam- 
ics that depend sensitively on the symmetry energy and 
its temperature and density dependence. 
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